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Abstract 

We compute the symplectic structure of the spin Calogero model in terms of 
algebro-geometric data on the associated spectral curve. 
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1 Introduction 



In a recent paper |T| I.M. Krichever and D.H. Phong made an important progress in 
the understanding of the symplectic structure of integrable models. Namely they were 
able to relate the canonical symplectic structure of the model with a naturally defined 
symplectic form on a suitable fibered space constructed with algebro-geometric data. In 
their paper this construction is illustrated by a number of examples. The next natural 
example to be treated is the spin Calogero model, and we show in this paper that the 
same general construction applies as well. This model has been solved in 0, where it was 
found that as compared to the scalar case a number of interesting new features enrich the 
algebro-geometric analysis. 



2 The model. 
2.1 Definition. 

The spin Calogero model consists of particles on a line with internal degrees of freedom. 
The particles are described by their positions Xi and momenta p,, together with spin 
variables fij. The Hamiltonian is: 

The non vanishing Poisson brackets for these degrees of freedom are: 

{fij,fki} = Sjkfil — ^iifkj (2.2) 



We see that the Poisson bracket is a Kirillov bracket and therefore is degenerate. We 



shall choose an orbit such that the matrix of spin variables / is of rank /. 
This system admits a Lax pair formulation, with Lax operator: 

Lijit, z) = pi6ij + (1 - 6ij)fij^{xi - Xj, z) (2.3) 

where for the elliptic model the special function $ is defined in terms of Weierstrass 
elliptic functions by: 

$(x,.) = ^fc^^'^^^^ 

It is elliptic in z and pseudo-periodic in x. It has been shown in |^ that the model is 
integrable only when all fa are equal, e.g. fa = 2. 



2.2 The spectral curve. 

In the article ^ it is noticed that this model can be related to the matrix KP equation, 
and this leads to its solution by algebro-geometric methods. In this paper we shall be 
concerned with the Hamiltonian aspect of the solution. The main object of interest is the 
spectral curve: 

R{k, z) = det {2kl + L{t, z)) = (2.4) 
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This defines a Riemann surface F of genus g = Nl — l{l + l)/2 + 1 (see 0) which is 
time-independent due to the Lax equation. The equation |2]J is polynomial in k of degree 
N with coefficients elliptic functions of z, hence presents F as an A^-sheeted covering of 
an elliptic curve with 2g — 2 branch points. On such a curve one usually defines a line 
bundle by taking at each point P = {z, k) the eigenspace of L{t, z) for the eigenvalue /c, 
i.e. the solutions of: 

i2kl + L(t,z))C{t,P) = 

In particular at the branch points the eigenspaces for the two colliding eigenvalues are 
generically of dimension 1, and L{t,z) is not diagonalizable for such z. 
Here special attention is required for the points above z = 0. 

2.3 The vicinity of z = 

At z = the Lax matrix L has an essential singularity since ({z) = 1/z + 0{z^). Hence 
one can write: 

L{z) = Diag(e^(^^-^°))L(2)Diag-^(e^(^'=-^°)) 

where Xq is some arbitrary origin, and L is meromorphic in a vicinity of z = 0. The 
eigenvectors are of the form C = (Ci) with: 

where Ci{P) is a locally analytic eigenvector of L. Since L{z) = —\{f — 21) + 0(1) 
we see that Ci{z = 0) is an eigenvector of f — 21 and the corresponding eigenvalues 
are of the form k = {—l + X/2)/z, where A are the eigenvalues of /, and iV — / of 
them vanish^. The class of functions having essential singularities at some points of the 
form ip = exp(a/z™)p(z) (here m = 1) with p locally meromorphic, and meromorphic 
otherwise are called Baker functions. The solution of many integrable models by algebro- 
geometric methods essentially depends on the construction of appropriate Baker functions 
on the spectral curve. A Baker function has properties similar to ordinary meromorphic 
functions, e.g. has the same number of zeroes and poles (consider the sum of residues of 
the regular differential dip/ip) but theorems such as Riemann-Roch have to be modified, 
for example there exists non trivial Baker functions with arbitrary prescribed g poles on 
a surface of genus g while for a meromorphic function one needs to prescribe g + 1 poles 
in the generic situation. To construct such a Baker function ip let us remark that one can 
find a unique normalized abelian differential of second kind uj2 and differential of third 
kind ujsipk) depending on the given poles, and the unknown zeroes pk of ip, such that 
dip/ip = U2 + uJz{pk)- By definition this form integrates to on A-cycles. Imposing the 
same condition on 5-cycles provides a system of g equations for the g unknowns pk- For 
a construction using theta functions see [Q]. Let us also remark that the quotient of two 
Baker functions with the same type of singularities is a meromorphic function to which 
one can apply the usual Riemann-Roch theorem. 

More generally one can define Baker line bundles, by imposing Baker conditions on 
the components around singular points, and this is the case for the eigenvector bundle of 

^ Remark that this defines genericaUy difi'erent branches of F which intersect at the singular point 
{z — Q,k — oo). After one blow-up I different points appear, and an ordinary multiple point of order 
N — / remains. An other blow-up at this point leaves us with N different points above z = 0. When we 
speak of T we have in mind this desingularization, and we speak freely of the N points above z = Q. 
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L. This line bundle is of Cliern class —{g—1). In fact in a nowhere vanishing section 
C(t, P) is constructed with g — 1 poles. 

Since Lij is not a symmetric matrix it is also convenient to construct an other Baker 
line bundle by considering the adjoint eigenvector equation: 

C+{t,P){2kI + L{t,z)) = 

which is also of Chern class —{g — 1). Here the Baker behaviour occurs with opposite 
exponent: 

C+(P) = e^(^'-^o)C'+(P) 

Note that these Baker line bundles are embedded in the ambient space of dimension 
N, which provides a pairing: 

<c+(p),c(p') >=J2cHPmp') 

i 

In particular for two points P and P' above the same 2;, corresponding to different eigen- 
values, one has < C+(P),C(P') >= due to the eigenvector equations, hence at the 
branch points < C+(P),C(P) > vanishes for any regular sections C"*", C. Moreover the 
singular factors exp— (xj — Xq) cancel above z = so < C^{P),C{P) > extends to a 
meromorphic function on F. 

Finally let Pa be the N points above the same z, and let Ca = C{Pa) be the corre- 
sponding eigenvectors. Except at the branch points, they form a basis of ambient space, 
and any vector V can be decomposed as: 

a=l <^'^a^'^OL > 



2.4 Remarkable abelian forms. 

The Riemann surface F presented by equation possesses several remarkable abelian 



differentials. First the form dz is well defined on F and has no pole. As such it has 
2g — 2 zeroes which are the branch points (where k is the local analytic parameter and 
dz/dk vanishes). The form of main interest is the form kdz which has poles only above 
z = 0. Finally let us take two non vanishing sections C^, C of the above line bundles. 
Then < C'^(P),C(P) > is a meromorphic function without any singularity above 2; = 
vanishing at the same points as dz hence VL = dz/ < C^{P), C{P) > is an other analytic 
abelian form vanishing at the 2g — 2 poles 7^ and 7^ of C and C'^. It will play an 
important role later on. 



3 The canonical symplectic structure. 



Due to the degeneracy of the Poisson brackets ^]2| one has to be careful about the choice of 
the symplectic variety. In fact this Poisson bracket is a Kirillov bracket for the coadjoint 
action of the group GL{N) acting on the dual of the Lie algebra gl{N) identified with 
itself by the invariant form {A,B) Tt{AB). The orbits are generically characterized 
by the eigenvalues of the matrix / which are in the center of the Poisson bracket. Here 
we shall consider matrices / of rank / with / different non-vanishing eigenvalues. Such an 
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orbit is of the form {g^^Ag\g G GL{N)} with A = Diag(Ai, . . . , A;, 0, . . . , 0). The tangent 
space at the orbit of / is the set of matrices U = [f, X] for any X G gl{n). In a basis where 
/ is diagonal this equation reads Uij = (Aj — Xj)Xij, hence Uij vanishes when Aj = Xj 
but is otherwise arbitrary. So the dimension of the orbit is 2NI — — I. At a point / 
the symplectic form on two tangent vectors [/, X] and [/, Y] is given by the well-defined 
formula: 

^;,([/,x],[/,y]) = Tr(/[x,r]) 

in accordance with Kirillov's prescription. Explicitly in a basis where / is diagonal this 
can be written: 

a;x = E ^^^^ (3.6) 

where by definition on a tangent vector [/, dfij{U) = Uij. 
Of course the symplectic form of our model is 



E dpi A dxi + ujk (3.7) 



The Hamiltonian is not invariant under the above GL{N) but it is is preserved 
by special subgroups. First we have the discrete subgroup of permutation matrices, i.e. 
the Weyl group, which simply operates by permutation of the N indices i 0, and more 
importantly we have the group of diagonal matrices, i.e. the Cartan torus, which operates 
by: 

fij di fijdj 

This action preserves the Hamiltonian which only depends on fijfji and all higher Hamil- 
tonians Tr(L"). The action of this toral subgroup induces a fibering of the orbits into 
fibers of dimension — 1 since multiple of the identity leave / invariant. The moment 
associated to this action is the collection of diagonal elements fa, that is A^ — 1 non trivial 
moments since on the orbit the eigenvalues of / are fixed, hence so is Tr(/). We consider 
the reduced dynamical system obtained by first fixing the moments to a common value 
fa = 2, and then quotienting by the stabilizer of this moment which is the whole diagonal 
group. This is the integrable system as considered above. 

We can now count the number of degrees of freedom. We have 2N degrees for the 
Xi,Pi, plus 2NI — — I for the orbit, minus 2(A^ — 1) due to the Hamiltonian reduction 
which ends up to a phase space of dimension 2(A^Z — l{l + l)/2 + 1) = 2g. 

It is a remarkable fact that the spectral equation is dependent on g non trivial 
integrals of motion. In fact it has been shown in p[ that considering only the order of the 
singularity at z = and knowing that all fa = 2 and / is of rank /, the spectral equation 
depends on g + 1 — 1 parameters. Here the eigenvalues of / are in the center of the Poisson 
algebra, and must not be counted as dynamical variables. There are / — 1 independent 
non vanishing eigenvalues of / since Tr(/) has previously been fixed to 2A^. Thus we end 
up with exactly g action variables. 

Our task is to find g other dynamical variables which will be used to construct the 
angle variables. It is known that such algebraically integrable systems linearize on the 
Jacobian of the spectral curve, hence it is natural to use g points ■jk on T as complementary 
variables. We have found above g — 1 poles of a non-vanishing section of the eigenvector 



^These actions extend obviously to actions on L by similarity transformations 
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bundle. If one defines C{P) such that Ci(P) = 1 as in 0, the poles of C are given by 
the vanishing on F of the first minor of the matrix L + 2kI. This is the algebraic equation 
which relates the 7^ to the dynamical variables of the system. Note that in this minor 
the variables Xi and pi have disappeared. This corresponds to reduction by translational 
symmetry, which leaves a phase space of dimension 2{g — 1). 

In order to get the full phase space it is convenient to take the product of such a 
section with an appropriate Baker function yielding a section with g poles and a fixed 
zero. It suffices to choose it with given zeroes canceling the poles of C, and one more 
zero at a fixed point Pq above z = 0. The singular behaviour above z = is taken to be 
expi^xo/z) so that we end up with Ci{P) proportional to exp(xj/z) having g dynamical 
poles, and similarly for C~^. Note that the poles 7^ of are determined when the poles 
7fc of C are given, since the abelian form Q = dz/ < C^{P), C{P) > is meromorphic with 
only singularity a double pole at Pq, hence has no residue at Pq, and one can choose the 
global normalization of such that: 

n = {\ + oii))dz p^Po 

z^ 

But such a form is uniquely determined when one fixes g zeroes 7^, and it has g other 
determined zeroes 7^ (otherwise the quotient of two such forms would be a meromorphic 
function with g poles and g other zeroes, which is generically forbidden). We summarize 
the definition of C and by stating their behaviour as 2; — *• 0: 

a(P) = e^Q(P), C+(P) = e^c+(P), z^^ (3.8) 

with Cj and regular (and vanish when P Pq). This together with the above condition 
on Q clearly determines C and up to a constant factor A on C and 1/A on when 
the 7fc are given, i.e. when the dynamical variables are given. Moreover, the q(P) for the 
points P above z = are N eigenvectors of / (for P = Pq of course we take Ci{P)/z 
as the corresponding eigenvector), and similarly for cf . 



4 The action— angle variables 
4.1 Some fiber bundles. 

Following the ideas of |I| we first introduce some natural bundles on the moduli space 
of our curves, or more specifically the ^f-dimensional space of action variables. Let us 
stress that we keep the eigenvalues of / fixed throughout our discussion. First we have 
the bundle Q whose fiber above a particular spectral function R{k, z) is the curve F of 
equation P = 0. Note that the differential dz is naturally defined on the fibers independent 
of R. We denote 6 the differential on Q, and we see that it can be splitted into a vertical 
part along the fiber Sy and an horizontal part which corresponds to differentiation with 
z fixed Sh- Hence if a^, i = 1, . . . , g are independent action variables we can write: 

. d . d . 

^ = 5^ * ^ 'dz 

In particular the two-form 6{kdz) = Snk A dz is well defined on Q. Note that k has 
simple poles on F whose residues are given by the eigenvalues of / as seen above, hence 
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the horizontal differentiation of these residues vanish, so the above form is regular on Q. 
This is one of the key observations in . 

Similarly, we define the fiber bundle whose fiber above R is the Jacobian of the 
curve of equation R = 0. A point of the Jacobian can be seen generically as an unordered 
set of g points 7^, i = 1, . . . , g of the curve F, and we also denote S the differential on J' 
which can be splitted as above, simply replacing the vertical part by: 

d 

where Zi is locally the 2;-coordinate of 7j. It follows that we can define [Q] on JT" a 
regular two-form u = S(J2i hdzi) which defines a symplectic structure. The notation 
is not accidental since we shall see that u precisely reduces to the canonical symplectic 
structure when we identify JT" to the phase space of our dynamical system. 

It is important to define, as in forms on J' with values abelian forms on the curve 
r associated to the corresponding base point of moduli space. In fact these forms are 
defined on the bundle whose fiber above a base point R is the Cartesian product of the 
curve r of equation R = and its Jacobian Jac(r). A generic point of such a bundle is 
described by g action variables Oj, a point P on F, and a set of g points 7^ of F. The 
differential 6 acts on the and the 7^, while the differential d acts vertically along F on 
P. 

We have found particularly convenient to introduce the meromorphic two-form on J': 

$ =< 6{C+{P)Q{P)) A 6L C{P) > (4.9) 

which due to the P dependence and the presence of fl{P) = dz/ < C^{P)C{P) > has 
values one-form on F. Here L is the Lax matrix depending implicitly on the and 7^, 
and the phase space differentiations appear in S{C~^Q) and 6L. Finally the brackets <> 
contract the i-indices of C, L and C~^. The main trick Q is to consider the sum of the 
residues of $ on F which must vanish. But {P)Q{P) is regular at 7^ for any value of 
the moduli and the 7^ hence so is its 5-differential. Hence the only possible singularities 
of $ are located at the g points 7^ of F where C{P) has a pole, and the points above 
z = where Q and 6L are singular. 

Writing these residues, we get a relation between a two-form on J' which happens 
to be the previously defined algebro-geometric symplectic form, and a two-form located 
above z = which boils down to the canonical symplectic form of our model. 

4.2 Algebro-geometric description of the canonical symplectic 
structure. 

We can now state the result of this paper: 



Proposition. 1 ) The canonical symplectic form |3.7| of the spin Calogero model can be 
written in terms of algebro-geometric variables as: 
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u = Y,SkiASzi (4.10) 

i=l 



A similar trick allows a considerable simplification of the proof of the main theorem in 
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where Zi = z{ji), h = kipa). 

2) If we take as angle variables the Abel transform of the divisor J^'^i namely the angles 
defined modulo periods ofT: 

dk=Y.r^k (4.11) 



where the ujk are a basis of regular abelian differentials dual to a basis {Ak} of A-cycles 
of r, then the canonically conjugated variables are given by: 

ai= i kdz (4.12) 

JAi 

Several remarks are in order at this point. First let us note that the first statement is 
exactly the same type of result advocated by Sklyanin in his famous solutions of various 
integrable models by his method of separation of variables |Q. In fact when one takes 
into account the above description of the poles of C one sees that it is really the same 
statement. However the result is obtained here without any appeal to the i?-matrix 
method. Moreover the second statement shows independently that the motion linearizes 
on the Jacobian of the spectral curve. 



Then let us note that formulae of type 4.12 have already appeared at various places. 



For example it is clear that they play a role in the classical and semi-classical analysis of 
the Neumann model 0, and have been introduced more generally in 0. In more recent 
analysis of quantum integrable models deformations of these formulae play a central 
role. One may hope that direct deformation of this description leads to the quantification 
of integrable systems in terms of algebro-geometric concepts. 

Finally, note that the proof of 1) =^ 2) is straightforward using a clever argument of 



In fact let us define ai according to equation 4.12, they are obviously moduli coordinates 



on the basis of our fiber space. Now as noted above 5{kdz) defined on Q is an analytic 
form on the fibers, hence can be expanded on the basis Uk of analytic one-forms. To find 
the coefficients, let us compute: 



d{kdz) daj ^ 
Ai (9a,- daj 



so we see that = tUk and we have: 



5{kdz) = ^ Sai A LJi 



Now, taking into account the value |4.11| of the angular variables we write the symplectic 
form as: 

uj = '^6{kj6zj) = ^6aiA uji{-fj) = ^ ^a^ A 69i 

j ij i 

which shows that the a, are indeed canonically conjugated to the 6',. One should note 
the simplicity of this derivation which essentially uses only one ingredient: the fact that 
the residues of the poles of kdz are killed by 6 so that 6{kdz) is regular. This can be 
contrasted to the involved computations inherent to previous approaches on this subject. 



7 



4.3 Proof of the main result. 



Let us compute the sum of the residues of the one-form $ on F (which are themselves 
two-forms on J\ As noted above we have to look at the points 7fc's and at the points 
above z = 0. In the vicinity of a pole 7^ we can write = [z — Z]-)"^ with \l/ regular 
for 2 = Zfc, hence: 

= + regular 

z — z-k 

(here bz^ = 6vZk), so the residue at 7^ is: 

T. ^ . <C+SLC>, , 
Res^,$ = -6vZk A ^ ^ (7fc) 

Note that 6h{C^Q) does not contribute to this residue because by definition SnZk = 0. 
But C+L = -2kC+ and LC = -2kC hence 6LC + L5C = -2k6C - 26kC which upon 
bracketing with yields < C^6L C >= —26k < C^C > so the final expression for the 
sum of the residues at the 7^ is: 

Res^j. = — 2 ^ 6kj A 6zj = —2uj 

3 

where lo is the previously defined symplectic form on J . Writing that the sum of the 
residues vanishes we have: 

2cu = 5^Resp,$ (4.13) 

a 

where the P^'s are the points above 2; = 0. 

Our next task is to compute this sum and show that it indeed reduces to the canonical 
symplectic form on our dynamical system. Let us remark that in the vicinity of z = we 
have a common local parameter z for all the curves in and is defined by the same 



normalization at Pq so is regular at Pq. Due to the normalization p.8| we can write: 

z 

We treat separately the contribution of the first term. The considered expression in the 
vicinity of P^ takes the form: 

dz < Cf5xi A SLijCj > . , 

~h~ < C+C > ^ ' 

Let us introduce the matrix Mij = 6xi A 6Lij which only depends on z{P). Summing on 
the N sheets we get, in view of equation the simple expression —dz/zTr{M), so the 



residue at z = is just J^i^Pi A 6xi. For the same reason the contribution proportional 
to 6Q can only pick diagonal terms in L and there is no residue. 

So we are left with the remaining contribution of e^6c^. Here we must examine the 
1/z terms present in SL. First note that all factors cancel between e^5c^, 5Lij and 
Cj. Then taking into account the identity: 

d 

— 0(x, z) = ^{x, z)[C{x + z) - C{x)] 



^ The point Pq presents no special problem since the factor present in the numerator cancels a 
similar factor in the denominator. 



8 



and the expansion 4>{x,z) = {—l/z + C{x) + 0{z)) exp{({z)x) one sees that all terms in 
6pi and Sxi are regular at 2; = 0, so only terms in Sfij contribute to the residue. Moreover 
(1 — Sij)6fij = 6fij because 6 fa = so one can simply replace 6Lij by {—l/z)6fij. The 
residue at z = is now obvious and the expression to be computed is therefore: 

<6c+A 6fco, > 
^ < r+r > 

where Cq is the value of c(P) at Pq,, i.e. they are the eigenvectors of /, and similarly for 
c+. 

To complete the calculation note that varying the equation /cq = AqCq, keeping A 
constant (recall that eigenvalues of / are fixed) one gets by contracting with the 
relation < c^Sfca >= (Aq, — A/j) < c^6ca >= — (A^ — A/j) < Sc^Ca >■ Expanding the 



above residue on the basis in view of one gets: 



^ < 6c+cp > ^ < c+6fcp > 



in which the sum can be restricted to Aq, 7^ A^ by the previous relation. But using ^ 
this is nothing more than: 



in which we recognize the expression of the Kirillov form uk- Finally we have proven 
that: 



2a; = ^ 6pi A 6xi + ujk 



5 Conclusion. 

The example treated in this paper shows once more the power and the generality of the 
method introduced by Krichever and Phong. Even in an intricate dynamical situation, 
the canonical symplectic form can be expressed in terms of algebro-geometric data with 
the same formula as in standard cases. Due to its simplicity and generality, this formula 
may very well extend to the quantum domain, as simple examples already indicate. 
Acknowledgements. We thank I. Krichever and F. Smirnov for useful discussions. 
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